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Abstract. The metallicity in portions of high-redshift galaxies has been successfully measured thanks to the gas observed 
in absorption in the spectra of quasars, in the Damped Lyman-o- systems (DLAs). Surprisingly, the global mean metallicity 
derived from DLAs is about l/10 th solar at 0SzS4 leading to the so-called "missing-metals problem". In this paper, we present 
high-resolution observations of a sub-DLA system at z abs =0.716 with super-solar metallicity toward SDSS J 1323-0021. This is 
the highest metallicity intervening quasar absorber currently known, and is only the second super-solar absorber known to date. 
We provide a detailed study of this unique object from VLT/UVES spectroscopy. We derive [Zn/H]=+0.61, [Fe/H]=-0.51, 
[Cr/H]=< -0.53, [Mn/H] = -0.37, and [Ti/H] = -0.61. Observations and photoionisation models using the CLOUDY software 
confirm that the gas in this sub-DLA is predominantly neutral and that the abundance pattern is probably significantly different 
from a Solar pattern. Fe/Zn and Ti/Zn vary among the main velocity components by factors of ~ 3 and ~ 35, respectively, 
indicating non-uniform dust depletion. Mn/Fe is super-solar in almost all components, and varies by a factor of ~ 3 among the 
dominant components. It would be interesting to observe more sub-DLA systems and determine whether they might contribute 
Q ■ significantly toward the cosmic budget of metals. 
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1 . Introduction ative to Zn probe the dust content and the relative abundances 

can also yield information about the nucleosynthetic processes 
Damped Lyman-a systems (DLAs) seen in absorption in the (e g ^ p ettini et al 1997; Kulkarni> Fall> & Trunffl) 1997. P6roux 
^ ; spectra of background quasars are selected over all redshifts in- et ^ 2002; et al 2004) A study of the cosmologica i 
dependent of the intrinsic luminosities of the underlying galax- evolution f the H 1 column density- weighted mean metallic- 
ies. They have hydrogen column densities, log N m > 20.3 and ity in DLAs (e g j Kulkarni & Fall, 2002) shows surprising re- 
are the major contributors to the neutral gas in the Universe at suks Contrary to most models of cosmic che mical evolution 
high redshifts (Storrie-Lombardi & Wolfe 2000; Peroux et al. (e g ^ Malaney & chaboyer, 1996; Pei, Fall & Hauser, 1999), 
2003b). But it has been suggested that at least some of the H 1 recent observations indicate at most a ^ evolution in DLA 
lies in systems with H 1 column density below that required by global met allicity with redshift for 0£z£ 4 (Prochaska et al. 
the traditional DLA definition, in the "sub-Damped Lyman-o< 2003; et al 2004; Ku]karni , et aL , 2 005; and references 
Systems (sub-DLAs)" with 19.0 < logAfoi < 20.3. The DLAs therein). Even theoretical models such as Smoothed-particle- 
and sub-DLAs offer direct probes of element abundances over hydrodynamics simulations (Nagamine, Springel, & Hernquist, 
> 90% of the age of the Universe. Zn is a good probe of the 2004) predict that the ^ DLA me tallicities could be 1/3 solar 
total (gas and solid phase) metallicity, in DLAs because Zn at z=2 5 and higher at lower redshifts. 
tracks Fe in most Galactic stars with [Fe/H]> -3, it is unde- 

pleted on interstellar dust grains, and the lines of the dominant Even at z=2.5, making a census of the predicted and ob- 
ionisation species Zn II are often unsaturated (e.g., Pettini et al. served neutral comoving densities of gas, Q, and metals, fi z , 
1999). Abundances of depleted elements such as Cr or Fe rel- one finds that most of the baryons are in the Lyman-ff forest 
but its metal content is extremely low. The measured value of 
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* Based on the UVES observations collected during the DDT ESO global mean metallicity is about 1/1 th solar. The metallicity 
prog. ID No. 274.A-5030 at the VLT/Kueyen telescope, Paranal, Chile of Lyman break galaxies is still poorly constrained; but, in 
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any case, these objects are known to be star-forming galaxies 
and may not be representative of the normal galaxy popula- 
tion. In total, these three components account for no more than 
a 10— 15% of what we expect to have been produced by z=2.5 
(Pettini et al. 2003; Bouche et al. 2005). The missing metals 
problem in low-redshift DLAs is even more surprising since the 
high global star formation rate estimates at z>1.5 (e.g. Madau 
et al. 1998) imply that higher metallicities should be expected 
at low redshift. 

It is possible that Qz(DLA) has been underestimated. 
Metal-rich DLAs could obscure quasars due to their possible 
high dust content (Fall & Pei 1993). This may be the reason 
for the apparent low metallicity in the DLAs observed in op- 
tically selected quasars (e.g., Fall & Pei 1993; Boisse et al. 
1998). Recently Vladilo and Peroux (2005) have shown that 
the fraction of high-redshift DLAs missed due to dust obscura- 
tion, could be up to 50%, which is consistent with the results 
of surveys of radio selected quasars (Ellison et al. 2001). They 
have estimated that at z~2.3, the real mean metallicity of DLAs 
could be 5 to 6 times higher than what is observed, which may 
help alleviate the missing metals problems. Indeed, systems at 
lower redshift may have significantly more dust at any given 
metallicity simply because the dust in these objects has had 
more time to process the elements. 

On the other hand, new lines of evidence are pointing to- 
ward lower Nui quasar absorbers like Lyman Limit Systems 
(LLS) and sub-DLAs being more metal-rich (Peroux et al. 
2003a; Jenkins et al. 2005). The dust bias, if real, is also likely 
to be less severe for metal-rich sub-DLAs as compared to the 
metal-rich DLAs due to the lower gas and therefore dust con- 
tent in the former, for a constant dust to gas ratio. Thus, the 
obscuration bias will affect the DLAs at a lower dust-to-gas ra- 
tio as compared to the sub-DLAs. This scenario is consistent 
with the recent radio surveys (e.g., Vladilo & Peroux 2005), 
but still needs to be further quantified observationally. Indeed, 
Zn measurements exist for only two sub-DLAs at low v 
the marginally super-solar sub-DLA toward Q0058+019 with 
z abs =0.61, log JVhi =20.08, and [Zn/H]=+0.08 (Pettini et al. 
2000), and the supersolar sub-DLA toward SDSS J1323-0021 
with z abs =0.72, logjV H i=20.21, and [Zn/H]=+0.40 (Khare et al. 
2004). Although Khare et al. reported [Zn/H]=+0.40 for this 
latter absorber, the modest resolution of their MMT data could 
not resolve the Mg I+Zn II A 2026 and Cr Il+Zn II A 2062 
blends. The extent of line saturation on the derived column den- 
sities was also unclear. With the goal of addressing these is- 
sues with high-resolution data, we obtained VLT/UVES spec- 
tra of this quasar, which are presented here. These new data 
are essential to confidently determine the metallicity by min- 
imising the problem of line saturation. Section 2 presents the 
observational set-up and data reduction process, while section 
3 presents the analysis. Section 4 provides a discussion of the 
results. 

2. Observations and Data Reduction 

Spectra of SDSS J1323-0021 (z em =1.390; SDSS mag 
g = 18.49) were acquired in service mode as Director's 
Discretionary Time (DDT) on 3 of March and 13 th of March 



2005 with the high-resolution UVES spectrograph mounted on 
Kueyen Unit 2 VLT (D'Odorico et al. 2000). Three exposures 
of length 4100 sec, 3500 sec and 4100 sec were obtained with 
standard 390+562 settings thus providing a wavelength cover- 
age of -3300A-4400A, 4700A-5600A and 5800A-6600A. 

The data were reduced using the most recent version of 
the UVES pipeline to accommodate for the new format of the 
raw fits file (version: uves/2. 1 .0 flmidas/1 .1.0). Master bias and 
flat images were constructed using calibration frames taken the 
closest in time to the science frames. The science frames were 
extracted with the "optimal" option. The spectrum was then 
corrected to vacuum heliocentric reference. The resulting spec- 
tra were combined weighting each spectrum with its signal-to- 
noise. The final spectra have resolution of 4.7 km s _1 at Zn II 
,12026. The spectra were divided into 100 A regions, and each 
region normalised using cubic spline functions of orders 1 to 5. 

3. Analysis 

Several lines of Zn II, Cr II, Fe II, Mn II, and Ti II were detected 
at z a b s =0.716. Fe I, Zn I, and Co II were not detected. The 
column densities were estimated by fitting multi-component 
Voigt profiles to the observed absorption lines using the pro- 
gram FITS6P (Welty, Hobbs, & York 1991) that evolved from 
the code used by Vidal-Madjar et al. (1977). FITS6P mini- 
mizes the^- 2 between the data and the theoretical Voigt profiles 
convolved with the instrumental profile. The atomic data were 
adopted from Morton (2003). 

The absorption profiles show a complex velocity struc- 
ture with a total of 16 components needed. The velocity and 
Doppler b parameters of the various components were esti- 
mated from the Mg I, Mg II, and Fe II lines. The component 
at 12 km s is negligible in most species except Mg I, Zn II 
and Cr II. The component at 194 km s is negligible in all 
species but Mg II. The component at -121 km s _1 is detected 
in Mg II, Fe II, and Mn II, but not in Mg I, Zn II, Cr II, or Ti II. 
The best-fit column densities in the individual components and 
their uncertainties were estimated assuming the same fixed b 
and v values for all species in a given component (Figure 0. 
The results of the profile fitting analysis are summarised in 
Table [0 Column densities in the few weak components that 
could not be well-constrained due to noise are marked with 
their contributions to the total column densities (listed in 
Table |5J are negligible. As an additional check, we also es- 
timated the total column densities using the apparent optical 
depth (AOD; Savage & Sembach 1991) method for the various 
detected lines and obtained results consistent with those from 
the profile-fitting method. The agreement was within 0.01 dex 
for Mg I and Mn II, within 0.1 dex for Fe II and Zn II, and 
within 0. 15 dex for Ti II. 

The Mg I /12026.5 contribution to the Zn II i2026.1 line 
was estimated using the component parameters for Mg I de- 
rived from the Mg I ,12852 profile. This contribution, indicated 
by a dashed blue curve in the top left panel of Fig^ was a small 
fraction of the observed strength of the /12026 line. The re- 
maining part of the ,12026 line was fitted with a 15 -component 
model for Zn II, using the same b values and velocities, but 
varying the column densities in the individual components. The 
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Table 1. Parameter fit to the z a b s = 0.716 sub-DLA model. Velocities and b are in km s 1 and N's are in cm 





Vel 


b 


Mgi 


Mgn 


Fen 


Znn 


Cm 


Mnn 


Tin 


[Fe/Zn] 


[Mn/Fe] 


N(X) 


-120.8 


7.1 




5.63el2 


5.82el2 






3.75ell 






0.79 


(T 








1.66ell 


4.08ell 






1.26ell 








N(X) 


-96.9 


16.3 




3.65el2 


3.64el2 


2.32el2 


<2.81el2 


2.97ell 


2.09ell 


-2.65 


0.89 


<r 








1.21ell 


4.31ell 


6.53ell 




1.93ell 


1.27ell 






N(X) 


-80.1 


6.2 


3.16el0 


3.55el2 


3.92el2 




<1.69el2 


2.09ell 


1.86ell 




0.71 


(T 






1.15el0 


1.36ell 


3.69ell 






1.23ell 


8.41el0 






N(X) 


-62.9 


5.5 


1.20ell 


8.13el2 


4.25el2 


<3.10ell 


<4.77ell 


2.01ell 


1.07ell 


-1.71 


0.65 


(T 






1.23el0 


3.68ell 


4.22ell 






l.Olell 


8.00el0 






N(X) 


-51.4 


5.0 


<1.33elO 


>2.86el2 


5.74el2 


<2.56ell 




1.44ell 


<6.92el0 


-1.49 


0.38 


<r 










1.30el2 






1.19ell 








N(X) 


-43.0 


7.4 


4.80ell 


>6.01el3 


2.41el3 


1.07el2 




1.98ell 


1.32ell 


-1.49 


-0.11 


(T 






2.29el0 




3.54el2 


5.18ell 




1.43ell 


9.85el0 






N(X) 


-13.0 


12.5 


3.04el2 


>2.38el4 


1.97el4 


2.57el2 




2.40el2 


4.79ell 


-0.96 


0.07 


(T 






1.48ell 




2.54el3 


6.22ell 




3.14ell 


2.16ell 






N(X) 


3.9 


30.5 


2.36el2 


>1.48el4 


4.00el4 






5.33el2 


<4.17ell 




0.10 


(T 






2.03ell 




2.63el3 






5.28ell 








N(X) 


12.0 


5.0 


1.44ell 






<3.66ell 


<3.20ell 




<3.89el0 




0.01 


(T 






4.81el0 


















N(X) 


28.3 


12.0 


2.06el2 


>7.14el3 


1.15el4 


1.53el2 


<2.37el2 


1.28el2 


3.80ell 


-0.97 


0.03 


(T 






1.25ell 




2.94el3 


7.60ell 




2.86ell 


1.53ell 






N(X) 


44.5 


11.0 


4.76el2 


>2.30el3 


2.89el4 


7.76el2 




4.34el2 


4.37ell 


-1.27 


0.16 


<r 






4.37ell 




5.42el3 


1.48el2 




2.90ell 


1.22ell 






N(X) 


56.6 


6.9 


8.70ell 


>4.62el3 


2.83el3 


<2.73ell 




6.30ell 


2.69ell 


-0.83 


0.33 


(T 






1.15ell 




1.82el3 






1.97ell 


9.48el0 






N(X) 


71.8 


11.9 


1.93el2 


>5.09el4 


3.12el4 


2.49el2 


6.07el2 


3.62el2 


1.58ell 


-0.75 


0.04 


cr 






5.19el0 




4.29el3 


7.30ell 


3.24el2 


2.22ell 


1.05ell 






N(X) 


92.7 


6.4 


2.42el2 


>2.99el4 


3.76el3 


7.35el2 


<1.54el2 


1.20el2 


2.09ell 


-2.14 


0.48 


cr 






8.51el0 




4.70el2 


1.49el2 




1.40ell 


7.87el0 






N(X) 


119.3 


8.0 


7.45el0 


4.49el2 


2.33el2 


5. Bell 


6.01el2 


1.98ell 


1.02e08 


-2.19 


0.91 


cr 






1.32el0 


4.90e07 


2.85ell 


4.48ell 


3.11el2 


1.27ell 


1.42e07 






N(X) 


194.1 


4.6 


2.44el0 


1.91el2 
















cr 






l.OOelO 


7.56el0 

















Zn II fit thus obtained was used to estimate the Zn II contribu- 
tion to the /12062 line. The remaining part of the /12062 line was 
fitted with a 15-component model for Cr II, using the same set 
of b values and velocities. This was the only available estimate 
for Cr II, since the Cr II AA 2056, 2066 lines lie in noisy regions 
and are undetected. The relatively large errors in the Cr II col- 
umn densities arise from the noisy nature of the /12062 line. 
The Fe II column densities in the components at low veloc- 
ities were constrained by using the Fe II /L12260, 2374 lines, 
since the stronger Fe II lines are saturated. The weaker Fe II 
components at high positive and negative velocities were con- 
strained in column density using the AA2374, 2382, 2600 lines, 
since these components are poorly constrained by the weaker 
X12260, 2374 lines. The Mg II ,12796, 2803 profiles were fitted 
together, but provide only a lower limit to Mg n owing to satu- 
ration in the central components. The relative abundances were 
calculated using solar abundances from Asplund et al. (2005), 
adopting the mean of photospheric and meteoritic values for 
Mg, Ti, Cr, Fe, Zn, and the meteoritic value for Mn. 

One concern is that sub-DLAs may be partially ionised in 
H, artificially enhancing the ratio of Zn II to H I, for instance. 
To investigate the ionisation corrections, we used the CLOUDY 
software package (version 94.00, Ferland 1997) and computed 



photoionisation models assuming ionisation equilibrium and a 
solar abundance pattern. We thus obtained the theoretical col- 
umn density predictions for any ionisation state of all observed 
ions as a function of the ionisation parameter U. Our findings 
confirm the observations: from a comparison of the observed 
and theoretical Mg II/Mg I ratios, we deduced that the gas 
in this sub-DLA is predominantly neutral (log U < -5) and 
the overall abundance pattern is probably not solar. This latter 
point is also clear from the relative abundances listed for each 
component in Table ^ It should be pointed out that there are 
no third ionisation stage detected in the system under study. 
Nevertheless, we do have Ti II, which has the same ionisation 
potential as H I. The fact that Ti is not suppressed compared to 
Fe or Mn also implies that the gas is neutral with low ionisation 
parameter. 

4. Discussion and Conclusions 

Table|2lists the abundances, using log Nhi = 20.21^2^1 (Khare 
et al. 2004) that we derived from Voigt profile fitting of the 
damped Ly-a line in the publicly available HST/STIS spec- 
trum of SDSS J1323-0021 (program GO 9382; PI: Rao). Rao 
et al. (2005) obtained log N m = 20.54+°^ from the same data 
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Fig. 1. Velocity plots of some of the element lines detected at Zzbs = 0.716 toward SDSS J1323-0021. The black tick-marks 
indicate positions of the components and the red curves denote the best-fit profiles. In the two top left panels, the upper set of 
blue tickmarks indicate positions of the components of Mg I /12026 and Cr II /12062, respectively, whilst the blue curves show 
their contributions to the absorption profiles. Note that different ordinates are used on the sub-frames. 



Table 2. Summary of total abundances. 



Id 


log iV, tai 


A(X/N) G 


[X/H]* 


Mgi 


13.26 ±0.01 






Mgn 


> 15.15 


-4.47 


> -0.58 


Fen 


15.15 ±0.03 


-4.55 


-0.51 ±0.20 


Znn 


13.43 ± 0.05 


-7.40 


+0.61±0.20 


Cm 


< 13.33 


-6.37 


< -0.52 


Mnn 


13.31 ±0.02 


-6.53 


-0.37±0.20 


Tin 


12.49 ±0.11 


-7.11 


-0.61±0.22 



* The error bars on [X/H] include the errors in log N(X) and log Nm- 

set. We use the former value since it gives a smaller residual 
with respect to the data and therefore regard this absorber as 
a sub-DLA. For either Nm, the strength of the Zn II lines de- 
tected in our UVES spectrum implies a super-solar metallicity. 
Using the standard definition: [X/H] = log[N(X)/N(H)] DLA - 
log[N(X)/N(H)] Q , we find [Zn/H]=+0.61. 

In principle, if Mg I A 2852 were substantially saturated, the 
contribution of Mg I A 2026 could be higher than our best-fit 
estimate. However, based on our apparent optical depth mea- 
surements and profile-fitting results, we estimate that it would 
take ~8 times more total Mg i than the value derived from the 
/12852 line to contribute the entirely of the /12026 line. Such a 
high value of Mg i can be ruled out by the observed profile of 
the /12852 line. (Of course, such a scenario would also be in- 
consistent with the observed strength and profile of the /12062 
line.) To understand this issue in more detail, we estimated the 



maximum Mg i in the dominant components that would still 
give the shape of the Mg I ,12852 profile consistent with the 
observed profile within the noise level in the continuum. This 
maximum total Mg i = 2.7el3 is about 50% larger than the 
best-fit value listed in Table 2. Putting this maximum Mg I 
model in the /12026 line, the corresponding total Zn n needed 
to fit the remaining part of the A 2026 line would be 2.5el3, 
lower by < 10% from our best-fit value. Thus [Zn/H] is at least 
> +0.59, indicating that our result would not be affected much 
by saturation of Mg I A 2852. Finally, as an additional check, 
we also estimated the maximum contribution of Cr II 2062 by 
rebinning our spectrum by factors of 10 or 20, measuring the 
upper limit for Cr UA 2056 in the rebinned spectrum. We then 
spread this upper limit for N(Cr II) over the 2062 profile, using 
the velocity model derived from the combination of the lines. 
We assume the same Fe/Cr ratio in all components, taking the 
percentage of Cr II in each components relative to total N(CrII) 
summed over all components to be the same as the correspond- 
ing percentage of Fe II in that component. Fitting the remaining 
part of the ,12062 line with a 15 -component model of Zn II, we 
estimated Az„n > 2.05el3, i.e. [Zn/H] > 0.50. 

The abundances of Fe, Cr, Mn, and Ti lie in the range of 
-0.4 to -0.6 dex, and indicate that this absorber is not only 
metal-rich, but also very dusty. Using the model from Vladilo 
(2004), we find that 95% of the Fe is in dust phase and the total 
metallicity is even slightly higher than 0.6 dex. This sightline 
also shows substantial reddening compared to the SDSS quasar 
composite (A(g - i) = 0.47; Khare et al. 2004). Considering 
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only the better-determined components between -13 and 100 
km s _1 , Fe/Zn varies by a factor of ~ 3 and Ti/Zn varies by 
~35. [Mn/Fe] varies by ~ 3, but indicates a super-solar Mn 
abundance with respect to Fe in all components. The relative 
abundance of Mn with respect to Fe is not expected to exceed 
the solar value as Mn, unlike Fe, has an odd atomic number. 
However, [Mn/Fe] >0 is often seen in the Galactic interstellar 
gas due to the stronger depletion of Fe on to dust grains. 

To summarise, our high resolution VLT/UVES data have 
allowed us to alleviate the saturation issue in Zn II and Cr II 
lines and therefore unambiguously prove the super-solar metal- 
licity of the sub-DLA at z abs =0.716 toward SDSS J 1323 -0021. 
If the dust obscuration bias for DLAs is indeed significant, 
as proposed by Vladilo and Peroux (2005), sub-DLA systems 
such as the one reported here could be better probes of dusty re- 
gions with significant past star formation (e.g. Lauroesch et al. 
1996; York et al. 2006), as similar DLA systems will be missed 
due to dust obscuration. On the other hand, it is also, possible to 
envisage a scenario where the dust obscuration is not very sig- 
nificant in quasar absorbers (as is indicated by the rising spec- 
trum of gamma-ray burst afterglows). In this scenario it is pos- 
sible that the sub-DLA systems may indeed be more metal-rich 
as compared to DLAs as indicated by observations of Peroux 
et al. (2003a) and by the observations of super-solar metallicity 
in such systems presented here and in Pettini et al (2003). With 
the large-scale spectroscopic surveys of quasars currently un- 
derway (e.g. the SDSS, York et al. 2000, York et al. 2001), such 
metal-rich sub-DLA systems may be found in large numbers. If 
future observations indeed find such systems, sub-DLAs may 
contribute significantly to the overall global metallicity. 
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